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A common external forcing can cause a saddle-node bifurcation in an ensemble of identical Duffing oscillators 
by breaking the symmetry of the individual bistable (double-well) unit. The strength of the forcing determines 
the separation between the saddle and node, which in turn dictates different dynamical transitions depending on 
the distribution of the initial states of the oscillators. In particular, chimera-like states appear in the vicinity of the 
saddle-node bifurcation for which theoretical explanation is provided from the stability of slow-scale dynamics 
of the original system of equations. Further, as a consequence, it is shown that even a linear nearest neighbor 
coupling can lead to the manifestation of the chimera states in an ensemble of identical Duffing oscillators in 
the presence of the common external forcing. 
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Ensembles of coupled oscillators are veritable black boxes 
which have been widely employed to understand a plethora of 
collective/co-operative dynamics observed in diverse natural 
systems Identification of simultaneous existence of a 

group of oscillators in harmony while the rest in dissonance 
in an ensemble of coupled identical oscillators has resulted in 
the notion of “chimera” states lH. Since the first identifica¬ 
tion of the coexistence of coherence and incoherence by Ku- 
ramoto and Battogtokh 0, many exciting developments have 
been made, including a number of experimental demonstra¬ 
tions IMl- A flurry of recent investigations revealed nonlocal 
coupling as a basic criterion for the onset of chimera in phase 
reduced models under weak coupling limit ||4]|5l, and also in 
systems far from the weak coupling limit ||9l. Very recently, 
the above limitation on the existing criterion for the onset of 
chimera has been liberalized to include an ensemble of glob¬ 
ally coupled oscillators CMni- Nevertheless, the mechanism 
for the birth of chimera states in an ensemble of identical os¬ 
cillators has been shown as the emergence of multistability in 
both nonlocal and global coupling configurations ll4ll5ll9 HT^ . 
However, it is not clear whether the converse is true, that is 
whether chimera can occur in a given multistable system ei¬ 
ther with nonlocal or with global coupling or even with much 
simpler couplings such as linear nearest neighbor coupling. A 
general criterion for the onset of chimera states in such dy¬ 
namical systems is yet to be unraveled. 

A periodically driven oscillator is the main and histori¬ 
cally first studied model in the classical theory of synchro¬ 
nization, where a triode generator is synchronized by a weak 
external periodic signal m. In reality, there exists a vari¬ 
ety of instances where the onset of such a collective behav¬ 
ior is invoked by an external force GHD. Examples include, 
phase locking of electrically decoupled spin torque nano os¬ 
cillators (STNOs) by an applied external microwave magnetic 
field ifBIfT^. polariton condensates in semiconductor micro¬ 
cavities that interact with the reservoir ini, frequency con¬ 
trolled devices, synchronization of micromechanical oscilla¬ 
tors using light QSl, etc. Here, we first consider an ensem¬ 


ble of identical bistable systems driven by a common exter¬ 
nal force without establishing any explicit coupling between 
the systems and investigate the underlying dynamics of the 
ensemble (or equivalently, that of an identical system for an 
ensemble of different initial conditions) as a function of the 
strength of the forcing. Then we investigate the effect of cou¬ 
pling leading to the manifestation of chimera states. 

The collective steady states resulting from the basin of at¬ 
traction of the ensemble of identical oscillators (or an ensem¬ 
ble of initial conditions) provides a clue about the basic crite¬ 
rion for the onset of chimeras in the ensemble of such systems 
with coupling. In particular, we find that if the basin of attrac¬ 
tion of a bistable (multistable) system leads to coexisting at¬ 
tractors of distinctly different nature under a common forcing, 
the ensemble of such systems will immediately lead to the on¬ 
set of chimera near the saddle node bifurcation of the ensem¬ 
ble under coupling. Eor appropriate forcing amplitude, the 
basin of attraction of the bistable system has essentially two 
domains of attraction. Initial conditions from one of the do¬ 
mains (coherent domain) lead to the same (chronotaxic mi) 
attractor while that of the other domain (incoherent domain) 
lead to quasiperiodic/chaotic attractors. Thus the collective 
steady state dynamics, whose initial conditions are uniformly 
distributed in the two domains, emerging from the basin of 
attraction resembles the chimera states under the common ex¬ 
ternal forcing. In the following, we will show as a general 
criterion that if the basin of attraction of a given dynamical 
system displays coexisting coherent and incoherent domains 
of attraction under the common external forcing, then even a 
linear nearest neighbor coupling can lead to the manifestation 
of chimera states in the ensemble of identical systems even for 
a low coupling strength in the same parameter space. This ap¬ 
proach will be extremely useful, and provides a cost-effective 
way in laboratory experiments to confirm the emergence of 
chimeras with just a single oscillator before experimenting 
with an ensemble of oscillators. In addition, the ensemble 
of identical oscillators driven by a common forcing opens up 
the possibility of analytical treatment which we will carry out 
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FIG. 1: (Color online) Distribution of the ensemble (N — 500) of 
Duffing oscillators in a double-well potential along with their average 
phase and average frequency distribution in the insets, (a) /e = 0, 
(b) /e = 0.088, (c) fe = 0.135, (d) /e = 0.15, (e) /e = 0.16 and 
(f) /e = 0.19. See text for explanations. 

in the following. 

In this Rapid Communication, we consider a bistable sys¬ 
tem with double-well potential as an individual oscillator ex¬ 
hibiting reflection symmetry which undergoes a pitch-fork bi¬ 
furcation EOll as a function of the system parameter. The basin 
of attraction of the ensemble of such identical oscillators leads 
to two cluster states corresponding to the bistable attractors 
with reflection symmetry. The common external force then 
leads to a saddle-node bifurcation of the ensemble by break¬ 
ing the reflection symmetry of the individual oscillators. For 
suitable values of the amplitude of the external force near the 
saddle-node bifurcation, we And the co-existence of coher¬ 
ent and incoherent steady states among the oscillator ensem¬ 
ble mimicking chimera states. Using a slow-scale approxima¬ 
tion and the method of direct separation 113, we separate the 
slow-scale dynamics from the original system of equations. 
The stability of the steady state solutions of the slow-scale 
dynamics provides an appropriate theoretical explanation for 
the observed dynamical transitions in the simulation. 

In general, an oscillator with a double-well potential of the 
form V (x) = x^ — x^ will have two stable fixed points and 
one unstable fixed point exhibiting bistability. A well known 
paradigmatic oscillator with this type of potential is the clas¬ 
sical model of Duffing oscillator ED, the ensemble of which 
with a common force is represented by 

Xi + axi — ujgXi + /3a;f — / sin(a;f) 

= fe sin(a;ef), z = 1,2,..., iV (1) 

where a = 0.5, Wq = 1,/3 = 1, w = 1 and / = 0.33 are the 
system parameters for which the individual Duffing oscillator, 
in the absence of the external force {fe = 0), exhibits periodic 
oscillations with period T « —, We = 0.5 is the frequency of 
the external forcing. We have fixed the number of oscillators 
as N = 500 and the system parameters as above throughout 
the manuscript. The Duffing oscillator has been investigated 
in great detail for its dynamical behavior with and without an 


additional external forcing for its chaotic nature and resonance 
phenomenon ll2Tll . 

Initial conditions of the ensemble of identical oscillators 
are uniformly distributed in both the wells with equal prob¬ 
ability in our simulation. Without any additional external 
force {fe = 0) the ensemble of individual oscillators will be 
oscillating synchronously in either of the wells as depicted 
in Fig. [TJa). However, the oscillators in both the wells ex¬ 
hibit out-of-phase oscillations with each other while oscillat¬ 
ing with the same frequency w = 1 (see insets of Fig. [TJa)) 
due to the symmetry property of Eq. Q- With the addition 
of the common external force to the ensemble, the synchro¬ 
nized oscillators in one of the wells desynchronize to form 
self organized clusters with distinct frequencies (see insets of 
Fig. [^b)), while the oscillators in the other well remain syn¬ 
chronized as shown in Fig.[2b) for fe = 0.088. Increasing the 
amplitude of the external force to fe = 0.135, the self orga¬ 
nized clusters desynchronize completely (Fig.[2c)). Thus the 
splitting of the basin of attraction of an ensemble of identical 
oscillators, for appropriate values of external forcing ampli¬ 
tude, into coherent and incoherent domains under the influ¬ 
ence of a common external forcing leads to the existence of 
chimera-like states. 

Upon increasing fe further, we And that the desynchronized 
oscillators hop back and forth to the other well where the syn¬ 
chronized group resides, the snap shot of which is depicted 
in Fig. [TJd) for fe — 0.15 again representing the chimera¬ 
like state. It is evident from Figs. [Jc) and [TJd) that the 
desynchronized group remains confined to one of the wells 
for fe = 0.135 and gets distributed in both the wells for 
fe = 0.15, thereby distinguishing two distinct chimera-like 
states C-I and C-II, respectively. The phase and frequency 
of the desynchronized group are clearly distinct from that 
of the synchronized group in the former case (see insets of 
Fig. [TJc)), whereas they are distributed about the phase and 
frequency of the synchronized group in the latter case as de¬ 
picted in the insets of Figs.[TJd) and (e). Some of the mobi¬ 
lized desynchronous oscillators are trapped in the other well 
to form a large cluster of synchronized oscillators for further 
larger values of the amplitude of the external forcing. For 
fe = 0.16, increase in the size of synchronized group is cor¬ 
roborated from the degree of phase and frequency distribution 
as in Fig. [^e). Finally, the desynchronized oscillators are all 
attracted to the other well to form a single synchronized clus¬ 
ter for /e = 0.19 as illustrated in Fig.[2f). Thus the ensemble 
of bistable oscillators loses its bistability to become monos¬ 
table by switching the stability of one of the stable fixed points 
via the chimera-like states as a function of fe- Existence of 
analogous chimera states between two populations with both 
inter- and intra-global couplings, as employed by Abrams et 
al ll22]| . has been demonstrated in discrete chemical oscillators 
based on the photosensitive Belousov-Zhabotinsky (BZ) reac¬ 
tion by Tinsley et al and in mechanical oscillator networks 
by Martens et al ||2l. Interestingly, we have observed similar 
transitions even in the states emerging from the basin of at¬ 
traction of an ensemble of identical bistable oscillators driven 
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FIG. 2: (Color online) Two parameter (/ — /e) phase diagram de¬ 
marcating various dynamical regimes. See text for details. 
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FIG. 3: (Color online) Analytical stability diagrams, (a) Profile of the 
fixed points of the slow-scale dynamical equation, Eq. l[^, (b) Level 
curve of Fig.|^a) depicting saddle-node bifurcation of the ensemble 
of oscillators for F — 0.65, (c) Level curve of Fig.|^a) for = 
0.0 showing pitch-fork bifurcation of the individual oscillators, and 
(d) Level curve of Fig. |^a) for F^ = 0.2 illustrating saddle-node 
bifurcation of the ensemble of oscillators. FPI, FP2 and FP3 refer to 
the three fixed point solutions of Eq. 0- 


by a common force. 

To obtain a global perspective, we have depicted the f — fe 
two parameter phase diagram in Fig. The parameter space 
corresponding to the collective steady state dynamical behav¬ 
iors discussed in Figs.[TJa)-(f) are labeled as (a)-(f) in Fig.|^ 
The steady states exhibit synchronous oscillations in either of 
the wells depending on the distribution of the initial condi¬ 
tions in the region marked as ‘SS\ Clusters of synchronous 
oscillators are seen in one of the wells in the region indicated 
by ‘CL’. Two distinct chimera-like states are observed in the 
regions denoted by ‘C — F and ‘C — IF, respectively. Syn¬ 
chronized single-well oscillations are found in the parameter 
space ‘S’. 

To explore the reason behind the observed transitions, we 
start with the dynamics of Eq. Q without any forcing. For 
/ = /e = 0,Eq.([TJ has two stable fixed points and one unsta¬ 
ble fixed point. Depending upon the distribution of the initial 
states, the ensemble of oscillators is attracted towards either 
of the stable fixed points. Eor finite values of /, the oscilla¬ 
tors corresponding to the stable fixed points oscillate in oppo¬ 
site phase but with the same frequency, observed in Eigs.[TJa) 
and l^a), in the asymptotic {t ^ oo) limit because of the 
invariant condition (xj,/sin(wf)) —?► (—/ sin(wt + tt)). 
Now, using the multi-time scale perturbation theory Ha and 
introducing a small parameter e << 1 as (a,wo,/3,/e) —>■ 
e(a, wq, £/3, sfe), the solution for Eq. Q can be separated as 
Xi{t) = yi{ujt) + Zi{et), where yi{Lut) = —■^sin{Lut) is the 
solution of Eq. ([^l for /e = 0 by neglecting the higher har¬ 
monics, and z{et) is a slow time scale perturbation. Plugging 
Xi{t) = yi{u}t) + Zi{et) in Eq. (|^ and averaging over 0 to 
we obtain 

Zi+az.i, + {F'^-ujl)z.i,+Pz^ = Fe, ( 2 ) 

where (x) = ^ / q “ xdt = z, F"^ = and Fe = 


{fe sm{ujet)) is non-zero for non-integer values of break¬ 
ing the reflection (left-right) symmetry leading to saddle-node 
bifurcation around which the dynamical transitions occur. 
Equation Q has three fixed points (a saddle and two stable 

foci/nodes) for Fe < Fe^± = The stability and 

the nature of the above fixed points determines the observed 
dynamical transitions discussed in Eigs.[2andJ^ Profiles of 
the fixed points of Eq. is depicted in Eig. |3[a) as a func¬ 
tion of F and Fe. The stable fixed points are symmetric about 
the saddle for Fe = 0 (see Eig. [^a)), whereas upon increas¬ 
ing Fe one of the stable (negative) fixed points and the sad¬ 
dle converge towards each other to merge at Fe = Fe^+ via 
saddle-node bifurcation as depicted in Eig. [^a). Conversely, 
the stable (positive) fixed point and the saddle can also con¬ 
verge towards each other to merge at Fe = Fe-. This means 
if F is replaced by —F, then we observe similar dynamics 
but in the other well. The fixed point profile for F = 0.65 
is shown in Eig. |^b) as a function of Fe. It is evident from 
both Eigs.j^a) and (b) that the separation between the saddle 
and the stable (negative) fixed point decreases monotonically 
upon increasing Fe, while the other fixed point remains un¬ 
altered. The degree of separation and the distribution of the 
initial states of the ensemble of oscillators altogether deter¬ 
mine the nature of dynamical transitions as a function of the 
forcing. 

Cluster states are formed in a range of Fe for fixed F de¬ 
pending on the degree of repulsion exerted by the saddle on 
the ensemble of oscillators whose initial states are distributed 
in the vicinity of the saddle and the stable fixed point. Such 
cluster states are represented in Eig.[2b). Indeed cluster states 
are formed in both the wells even for Fe = 0 (fe = 0) but for 
large values of F (/) as can be seen in Eig.[^ However, states 
mimicking chimera do not exist for Fe = 0 because both the 
fixed points approach the saddle resulting in complete inco- 
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FIG. 4: (Color online) Spatiotemporal plots (left column) and prob¬ 
ability distribution of time averaged frequency (right column) of an 
ensemble of Duffing oscillators with linear nearest neighbor coupling 
in the presence of the common external forcing depicting (a)-(b) clus¬ 
ter states for e — 0.06 and = 0, (c)-(d) chimera-like states for 
£ = 0.0 and /e = 0.14, and (e)-(f) chimera states for e = 0.06 and 
/e = 0.14. 


herent oscillations in both the wells for large F ED. Finally, 
the fixed points collide with the saddle to create a new sta¬ 
ble fixed point through pitch-fork bifurcation (see Fig. [^c)) 
resulting in double-well oscillations. In contrast, for a finite 
value of Ff., above Fe — 0, the ensemble of oscillators un¬ 
dergo double-well oscillations via cluster, chimera-like states 
and single-well oscillations, where the symmetry between the 
fixed points are broken and saddle-node bifurcation occurs at 
the onset of single-well oscillations as a function of F as de¬ 
picted in Fig.|^d) for Fe = 0.2. 

As the separation between the saddle and the stable fixed 
point narrows down further on increasing Fe beyond that of 
the cluster states, the saddle expels the oscillators away from 
the stable fixed point resulting in incoherence among the os¬ 
cillators in the well, while the oscillators in the other well are 
in coherence by which C-I emerges in Fig. [^c). For further 
large values of Fe the separation between the saddle and the 
stable fixed point becomes negligibly small and they eventu¬ 
ally merge together to form a saddle-node bifurcation. In this 
narrow range of Fe (as seen in the two parameter phase dia¬ 
gram in Fig. 1^, some of the oscillators whose initial condi¬ 
tions lie in the vicinity of the saddle makes round trip among 
both the wells and some of them are attracted in their opposite 
well depending on the value of Fe elucidating the existence 
of C-II states, which we have observed in Figs[2d)-(e). Be¬ 
yond the saddle-node bifurcation, the saddle and stable fixed 
point disappear only with the existence of the stable positive 
fixed point in the other well to which all of the oscillators are 
attracted finally for sufficiently large Fe confirming our dis¬ 
cussion on Fig.[2f). 

Now, we will show that even a linear nearest neighbor cou¬ 
pling can lead to the manifestation of chimera states in an en¬ 
semble of Duffing oscillators in the presence of the common 


external forcing represented as 

Xi + axi — ujqXi + Px^ — f sin(a;f) 

+ eG = feSm{oJet),i = 1,2,...,N 

(3) 

where £ is the coupling strength and G is the coupling 
function characterizing the linear nearest neighbor coupling 
{xi+i + Xi-i — 2xi) while the system parameters are the 
same as in Eq. Q. Spatiotemporal plots and the probabil¬ 
ity distribution of the time averaged frequency (17) are de¬ 
picted in the left and the right columns of Fig. for various 
values of the system parameters. In the absence of the com¬ 
mon external force fe = 0 the ensemble of Duffing oscilla¬ 
tors with nearest neighbor coupling displays cluster states as 
shown in Figs.Qa) and (b) for the coupling strength e = 0.06. 
The ensemble will exhibit chimera-like states, as discussed 
above, when driven by the common external forcing in the ab¬ 
sence of any couplings {e = 0.0) between the oscillators (see 
Figs. Qc) and (d)) for fe = 0.14. Chimera states manifest 
even due to a linear nearest neighbor coupling in the pres¬ 
ence of the common external forcing as shown in Figs. Qe) 
and (h). The spatiotemporal plot and the probability distri¬ 
bution of (17) is depicted in Figs. |^e) and (f), respectively, 
for the coupling strength e = 0.06 and the forcing amplitude 
fe = 0.14. It is clear from these figures that there is a si¬ 
multaneous emergence of synchronized and desynchronized 
domains among the ensemble of identical Duffing oscillators 
elucidating the existence of chimera states due to the near¬ 
est neighbor coupling in the presence of the common external 
force in the same parameter regime where chimera-like states 
are observed without coupling. The effect of nearest neighbor 
coupling is clearly evident from Fig.j^f), where there is a shift 
in the frequency of the synchronized domain compared to that 
of the chimera-like states in Fig. |^d), while the frequency 
of the desynchronized domain has a much wider distribution 
than that of the chimera-like states. Thus it is evident that the 
chimera state manifests due to the nearest neighbor coupling 
but in the presence of the common external force. Chimera 
states emerge even for such a low value of coupling strength 
£ = 0.06, essentially because of the precursor chimera-like 
states in the absence of coupling. The above results indicate 
that even an ensemble of oscillators with linear nearest neigh¬ 
bor coupling can exhibit chimeras under suitable conditions. 
Note that the system can also be interpreted as an ensem¬ 
ble of Duffing oscillators with double forcing and with lin¬ 
ear nearest neighbor coupling. We have further verified the 
generic nature of our results in ensembles of Stuart-Landau 
oscillators and Rossler oscillators in the presence of the com¬ 
mon external forcing, the details of which will be published 
elsewhere. 

To summarize, surprisingly we have found that the basin 
of attraction of an ensemble of identical bistable oscillators 
can indeed display coexisting coherent and incoherent do¬ 
mains with distinctly different nature of attractors mimicking 
chimera states under a common forcing, even without any ex¬ 
plicit coupling, using Duffing oscillator as a typical example. 
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The common force facilitates a global saddle-node bifurca¬ 
tion of the ensemble by breaking the symmetry of the individ¬ 
ual oscillators exhibiting pitch-fork bifurcation in the absence 
of the external forcing. The spontaneous splitting of the en¬ 
semble is found to emerge near the saddle-node bifurcation 
under an appropriate forcing. We have also provided an ap¬ 
propriate analytical treatment of the original system of equa¬ 
tions, where we have separated the slow-scale dynamics and 
investigated the stability of the underlying fixed points, which 
provides necessary explanation confirming the observed dy¬ 
namical transitions in the simulation. Further, it is also shown 
that even a linear nearest neighbor coupling can lead to the 
manifestation of the chimera states in an ensemble of identical 
Duffing oscillators but in the presence of the common external 
forcing. As driving/influencing others is a natural behavioral 
tendency in ecology such as in a herd of sheep, flock of birds, 
colony of ants, hive of bees and many more, as well as in epi¬ 
demics, in neuroscience, in social networks, etc., there lies 
every possibility that such an emergent behavior may exist in 
several natural systems. Our results may open up potential 
activities in the identification of chimera states in appropri¬ 
ate natural systems. In particular, our results have elucidated 
that still there are new avenues open to extend the horizon of 
the existence criteria for the framework of chimera with more 
simple and realistic couplings/situations that fits with the real 
world examples. More significantly, our results will serve as a 
basic framework to ensure the existence of chimeras in labora¬ 
tory systems before performing experiments using ensembles 
of such systems. 
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